Invading interfaces and blocking surfaces in high dimensional disordered systems 
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We study the high-dimensional properties of an invading front in a disordered medium with random 
pinning forces. We concentrate on interfaces described by bounded slope models belonging to the 
quenched KPZ universality class. We find a number of qualitative transitions in the behavior of the 
invasion process as dimensionality increases. In low dimensions d < 6 the system is characterized by 
two different roughness exponents, the roughness of individual avalanches and the overall interface 
roughness. We use the similarity of the dynamics of an avalanche with the dynamics of invasion 
percolation to show that above d = 6 avalanches become flat and the invasion is well described as 
an annealed process with correlated noise. In fact, for d > 5 the overall roughness is the same as 
the annealed roughness. In very large dimensions, strong fluctuations begin to dominate the size 
distribution of avalanches, and this phenomenon is studied on the Cayley tree, which serves as an 
infinite dimensional limit. We present numerical simulations in which we measured the values of the 
critical exponents of the depinning transition, both in finite dimensional lattices with d < 6 and on 
the Cayley tree, which support our qualitative predictions. We find that the critical exponents in 
d = Q are very close to their values on the Cayley tree, and we conjecture on this basis the existence 
of a further dimension, where mean field behavior is obtained. 



I. INTRODUCTION 

The problem of interface growth and wrinkling in dis- 
ordered and noisy systems has been the subject of much 
interest ||]. A special attention was given lately to in- 
terface growth in disordered systems. Examples of such 
systems are two fluid displacement flow in porous media 
PI , invasion of water into paper Isl and magnetic domain 
movement in spin systems with quenched disorder. The 
distinction between such disordered systems, where the 
noise is quenched, and cases were the motion is induced 
by random annealed noise has been an important part 
of this study (llMl. It is known for some time that the 
critical properties of interfaces in the two kinds of sys- 
tems are different |l|,||,||j|]. As a consequence of the 
time independence of disorder in a quenched system an 
interface invading a quenched medium may be pinned, 
if the driving force is smaller than some threshold value, 
while in the annealed version it will always grow. The de- 
pinning transition is a critical phenomenon, in which the 
correlation length diverges as one approaches the thresh- 
old driving force value. Thus, there is usually a marked 
distinction between the behavior of the quenched noise 
problem and the annealed versions, and they define dif- 
ferent universality classes. 

Advancing interfaces are usually characterized by a 
roughness exponent x relating the average width of the 
interface W to the system size L, 



W ^L^. 



(1) 



In most cases, the interface has also a self-affine struc- 
ture characterized by the same roughness exponents. 



A particular difference between annealed and quenched 
growth is that the roughness exponent Xq of an inter- 
face growing in quenched disorder is generally different 
from Xa, the roughness of the annealed counterpart. An 
important example is bounded slope growth, which may 
be described by the KPZ equation for the annealed 
case |l|. The quenched model pIpPH has a different 
roughness. 

An interesting aspect of quenched growth processes is 
the mapping which exists between such problems and the 
percolation of random surfaces. In fact, a pinned inter- 
face describes a blocking surface which traverses from 
one side of the system to the other. If the interface is 
depinned, on the other hand, no such blocking surface 
exists. In the simplest example of 1-1-1 dimensions, the 
blocking surface describes a line in a directed percola- 
tion cluster |§,||. In other dimensions such generalized 
percolation processes are less well understood. Thus, we 
can identify the directed percolation transition with the 
depinning transition in 2 dimensions. 

We want to emphasize that the analysis presented in 
this paper relates to a universality class which is realized 
in many different contexts. It was found in studying other 
systems, that much useful information and insight can 
be gained by studying the behavior in spaces of varying 
dimensionality. This is the case also for the interface 
growth model studied here. Since the model is rather 
involved we discover a number of qualitative transitions 
as the dimensionality increases. 

Turning back to the dynamics of interfaces we note 
that critical behavior in such interfaces can be observed 
using two different driving mechanisms. First, one ob- 
serves critical properties near the threshold force level 



Pc, required for the interface depinning. Alternatively, 
one may drive the system at the point of weakest resis- 
tance using a constant infinitesimal current. This driving 
causes the system to self-organize into a critical state. 
The critical interface in both cases is that of a self-afRne 
interface with infinite range correlations. Interfaces pro- 
duced by the two methods are intimately related. We 
will use both methods in our discussions. 

Analyzing interface growth driven by infinitesimal cur- 
rent, we can define a key concept for the understanding 
of such processes, namely the associated process (APs) 
|^,||. Associated processes are defined as the domains 
covered by the growth at periods during which the resis- 
tance of the medium is below a certain threshold value. 
In annealed disorder one cannot define such processes 
and this is the main source for the distinction between 
annealed and quenched growth models. We analyze the 
statistical properties of the APs, and in particular their 
accompanying roughness exponent Xc which is defined by 
the scaling of the height of the APs with respect to their 
lateral dimension. We stress that Xc is an independent 
exponent and may be different from Xqi since the latter 
is created by a process of stacking of independent critical 
APs. In fact, Xc is only a lower bound for the overall 
roughness Xq [1|- Thus in general, blocking surfaces have 
two independent roughness exponents. The first is the 
global one, and the second is the scaling exponents of 
the voids between the surfaces. 

In a previous letter [nO[ we extended the discussion 
and analysis of interface growth as a series of APs to an 
arbitrary number of space dimensions, and used the AP 
concept to derive a scaling theory of interface growth 
in all dimensions. We argued from general principles 
that the APs become less rough as the dimensionality 
increases and eventually become flat objects at a finite 
dimension dc. For space dimensions d > dc the APs are 
fractal. Note that this dimension is critical in the sense 
that the scaling relation change qualitatively reflecting a 
fundamental change in the dynamical process. However 
a second critical dimension can exist where the critical 
behavior will be the same as in the infinite dimensional 
limit. 

In this paper we augment the discussion in |Ig] in sev- 
eral ways. We consider the dynamical process that oc- 
curs during the evolution of a critical AP. The set of sites 
which are active at a certain instant is sparse. Therefore 
we conjecture that it is a percolative dynamical process. 
Using this claim we can put limits on the maximal pos- 
sible AP roughness Xc as a function of d and we are able 
to show that dc = 6 for the quenched KPZ universality 
class. This is related to the observation in m^ that the 
dynamical exponent z equals 2 for d > 6. 

We next present simulations of the model on lattices of 
up to six dimensions, verifying that the dc is indeed equal 
to 6. We analyze the the overall interface growth and 
show that the interface roughness exponent, Xq are the 



same as the exponents Xa of the annealed KPZ equation 
above d = 4. 

The infinite dimensional limit of this problem can be 
analyzed using a realization of the model on a Cayley 
tree with an additional height coordinate. Section three 
is devoted to this analysis. We performed simulations 
on this model and found two main outstanding features. 
The first one is the existence of strong anomalous behav- 
ior, resulting in power law distribution functions with 
changing exponents in the subcritical range. Secondly, 
near the depinning transition a critical behavior emerges 
with a strong anomalous background. The values of the 
scaling exponents seem to be the proper infinite dimen- 
sional limit of the model. In fact the scaling exponents in 
six dimension, which is the highest dimensional system 
we analyzed numerically, are quite close to the results on 
the Cayley tree. This is an indication that the critical di- 
mension of the problem, where all exponents reach their 
infinite dimensional limit is close to 6. 

In section 3 wc discuss two variants of the model on 
the Cayley tree. The first is a straightforward general- 
ization of the bounded slope model we just described, in 
which the bounded slope condition permits avalanches 
propagating in all directions, allowing for self-interaction 
by a re-activation of a site by a daughter site. This self- 
interaction prevents us from obtaining a closed form so- 
lution. Hence we define a simplified version in which the 
bounded slope condition is enforced only in a preferred 
direction. This directed version is simplified enough that 
analytic calculation can be done, using a mapping to an 
infinite state branching process. However, the directed 
version is not a satisfying high-dimensional limit, since 
the anomalous fluctuations mentioned above swamp out 
the critical behavior of the transition, and do not allow 
for an explicit calculation of scaling exponents. 



II. THE INVASION PROCESS AND INTERFACE 
STRUCTURE IN HIGH DIMENSIONS 

1. Basic definitions and scaling relations 



The discrete Buldyrev-Sneppen model [g|j5[|7|], which 
belongs to the quenched KPZ universality class, was orig- 
inally defined on a l-fl dimensional lattice. At each site 
there is a quenched random pinning force f{s) distributed 
uniformly between and 1. There are two related ways 
to define the invasion process. In the constant current 
version one has to locate the site of minimal / value on 
the interface. This site is activated by increasing the in- 
terface height at that point. Next one checks whether the 
activation of this minimal site creates a local slope which 
is larger than a threshold value. If such an event occurs, 
the neighboring site is also activated to reinstate a sub- 
threshold local slope. These two operations are carried 
out repeatedly creating an advancing interface. Another 



algorithm for advancing the interface is the activation of 
all sites on the interface whose value of pinning force / 
is smaller than a fixed driving force value /o, followed 
by avalanches to enforce the bounded slope condition, 
and repeating these operation. If /o is chosen such that 
it is smaller than a critical value Pc, the interface will 
be pinned eventually, and if /o > Pc, the interface will 
continue to propagate indefinitely. 

A very fruitful way of describing the invasion process 
by constant current is to introduce the concept of associ- 
ated processes l§,0]. They are defined by examining the 
value of / for each activated site. An /o AP is a set of 
activations with / values below the threshold value /q. 
An activation above /o concludes an /o AP and starts a 
new one. The APs of /2 may comprise several /i APs if 
/i < /2- The largest APs are the Pc APs. In the constant 
current algorithm, no site with f > Pc can be activated. 
An alternative way to construct an AP, is to switch to 
the constant force /o algorithm immediately after updat- 
ing a site with / > /o- The process will be blocked, and 
the cluster which is created is precisely an AP, and this 
fact relates the two driving mechanisms. 

The definitions of the model and the two driving meth- 
ods can be carried through unchanged to lattices of 
higher dimensionality. One defines an interface height 
coordinate in each site, and performs the invasion pro- 
cess according to the rules just describes. The lattice 
does not have to be an Euclidean lattice, and we will 
consider interfaces growing on a Cayley tree in section 3. 

For /o close to Pc the APs display critical behavior, 
where A/ = pc—fo measures the closeness to criticality of 
the AP. We are interested in the distribution of AP sizes 
(or masses) s defined as the number of activations in the 
AP, and in the RMS lateral dimension r|| . The following 
scaling forms are known to be valid for d-dimensional 
APs near fc'. 

Pis) = s-^gis/Af-n, Pin) = r^'"g||(r||/A/-''ii), 

(2) 

and the fractal dimension D of the APs is defined by 

s^r^. (3) 

The scaling relations 



1 



z/||(d+l-r||). 



(6) 
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1 = DiT ~1) 



(4) 



follow immediately from the definitions (||) and (^. If 
D > d, the APs are rough objects, that is, they have a 
width r±_ related to ru by 



r± 



II 



(5) 



where Xc = D — d. 

An exponent relation proved in H using the properties 
of the model is 



This relation reduces the number of independent static 
exponents to 2 and remains valid as long as the AP rough- 
ness exponent Xc is positive. A relation equivalent to (H) 
is 



7 = 1 + z^Xc 



(7) 



where 7 — 1^(2 — r), is the scaling exponent of the average 
cluster size with respect to A/. 

We will show below that there exists a finite dimension 
dc such that Xc = when d > dc- Above dc the scaling 
relation (0) becomes simply 7=1. However, the other 
exponents and the fractal dimension continue to depend 
on the dimension, and therefore dc does not mark the 
transition to mean field behavior. 

The geometric structure of an AP is qualitatively dif- 
ferent in d > dc. Studying the rules of the model in- 
dicates that whenever Xc is positive, the APs cannot 
be fractal, since rough clusters have a typical width rj_, 
which through avalanches prevents the formation of holes 
with diameter smaller than rj_. For d < dc the APs may 
be multiply connected, but not fractal. Above dc, since 
D < d the AP become objects with order 1 thickness, 
which may fractalize by having holes of any size. 



2. Three types of roughness exponents 

In the previous subsection we defined the roughness ex- 
ponent Xc which characterizes the statistics of the APs. 
It is important to realize that Xc is not necessarily equal 
to the exponent Xg defined in eq. (|lj), which describes 
the interface as a whole. However, a simple considera- 
tion shows that Xg cannot be smaller than the value of 
Xc- In fact, from time to time an AP is generated which 
encompasses the whole system, and then the entire in- 
terface identifies with a critical AP (see also g]). The 
interface growth is also related to the annealed process, 
since the stacking of APs is similar in spirit to the pro- 
cess described by annealed equations. This issue is dealt 
with in more detail below, but here we would like to note 
that it follows from our discussion that the roughness of 
the annealed dynamics Xa is also a lower bound for Xq: 
so we may write 



Xq > max(xa,Xc)- 



(8) 



Since the quenched growth includes a non-trivial interac- 
tion between the processes defining Xa and Xc its rough- 
ness can actually be larger than both x's. 

As an example one can consider the 1 -f 1 dimensional 
case. Since the APs are related to directed percolation, 
their roughness is Xc — 0.63 P,pi. The roughness for the 
annealed dynamics is Xa = 0.5, so that Xc dominates Xa- 
In fact, it has been found (see B) that in this case the 



system roughness Xq is slightly larger than Xc- In higher 
dimensional systems Xc decreases, and the difference be- 
tween Xq and Xc becomes more pronounced. 

As an application of this distinction, we note that when 
constant force driving is used, starting from a flat inter- 
face as initial conditions, the created clusters will have 
roughness Xq rather than Xc- In fact the clusters will 
consist of many APs since the initial surface cuts through 
them. This difference will be very manifest in high di- 
mensions, where the APs arc flat and fractal whereas the 
constant force clusters are rough pl]|. 



3. Dependence of Xc on dimensionality 

It is apparent that growth in APs can happen in 
two ways; either upward by activations, or sideways via 
avalanches. The sideways action becomes more impor- 
tant as the dimensionality increases, and Xc therefore 
decreases as a function of d. A naive expectation based 
on this scenario is that in sufSciently high dimension the 
sideways growth becomes non self-interacting. If this ar- 
gument were true, the APs should belong to the percola- 
tion universality class in large enough dimensions. This 
simple behavior does not occur, since the bounded slope 
condition induces self-interaction which does not decrease 
when d increases, as we explain in detail in the next sec- 
tion where we analyze the model on the Cayley tree. 

There is however a subtle connection between the dy- 
namics of a growing AP and the dynamics of a percola- 
tion cluster. In the following subsection we exploit this 
analogy to argue that in some finite dimension dc the APs 
become flat objects with Xc — 0. As explained above, 
dc does not mark the transition to mean field behavior. 
Some exponents, such as the AP fractal dimension _D, 
change as a function of d above dc. Hence, for d > dc the 
APs are fractal objects with D < d (see subsection 2.1). 
We expect that D saturates at some finite value Doo, as 
indicated by the Cayley tree analysis presented in the 
next section. 



4. Dynamical scaling and critical dimension 

The discussion presented so far has concentrated on 
the AP static exponents such as r and D. It is inter- 
esting that considerations involving dynamic exponents 
provide further information on static exponents as well. 
Although the arguments given here are not rigorous, we 
believe that they still give a consistent and informative 
picture of the invasion process. 

The dynamical exponent z is defined as follows: Con- 
sider an AP initiated at some site following its growth 
activating at each time step all the available sites at that 
moment. The lateral size of the AP grows as a power law 
in time 



t 



(9) 



It is easy to show that the fractal dimension of the set of 
activated sites at any instant is 



Dae = d + Xc 



(10) 



Consider the ongoing process of the activations on the 
surface of the cluster. Since the fractal dimension Dae is 
smaller than d, (since Xc < ^ < z) the set of active sites 
is so sparse that the time interval between subsequent 
activation of a single site will be very long (in fact we 
can show that it scales as r^X'=/'^^~^'>). This indicates that 
the dynamical process of the AP growth is a percolative 
process, similar to the process the dynamical invasion 
into a percolation cluster. The first conclusion from this 



argument is that z 



the dynamical exponent of 



invasion percolation. This relation has been suggested in 
|12| and verified numerically in dimensions d < 6. 

The AP dynamics in fact resembles a set of invasion 
percolation processes occuring simultaneously, so that 
Dae may be larger than the dimension of the set of active 
sites in a percolation process Dp^r — Zper- However Dae 
cannot be larger then the fractal dimension of a full per- 
colation cluster Dper without completely destroying the 
dynamical picture. Hence consistency demands that 



Dae = d + Xc - z < D.^ 



and thus, 



Xe Jl ^per i Zpcr 



per : 



d<6-d. 



(11) 



(12) 



In the last inequality we used the largest possible values 
for Dpcr and Zper which are the mean field values. The in- 
equality (|l^) gives an upper bound for the AP roughness 
Xe in each dimension. We verified that this bound in- 
deed holds for the numerical values of Xc, and is actually 
rather sharp (see table 1 in subsection 2.5). Moreover, 
(n2) implies that when d > 6, Xe must be 0, so we have 



de<6 



(13) 



We note that for d > dc our arguments become irrele- 
vant since the AP is then a fractal objects and its growth 
in this case will be occur mostly in the boundary and 
not in on its bulk. Therefore we expect that the pro- 
cess is no longer a percolative one. It is amusing to note 
that for this model only the low dimensional processes in 
d < dc — 6 are percolative. 



5. Numerical simulations in high dimensions 

We performed numerical simulations of the Buldyrev- 
Sneppen model and compared them with the previous 
theoretical scenario. The simulations were of either of 
two kinds: Interface growth by constant current on a 



hyper-cubic lattice with periodic boundary conditions, 
and invasion of a single cluster in a Cayley tree as an 
effective infinite dimensional lattice. 

In each of the finite-dimensional simulations the num- 
ber of activations was of the order of lO^*', and thus sta- 
tistical errors were rather small. However, in high dimen- 
sional systems we were severely constrained by memory 
considerations. For example in 6 dimensions, the highest 
dimensionality in our simulations, the hyper-cube mea- 
sured only 16 sites in each directions. Thus, finite size 
effects were very strong in the high dimensional cases, 
leaving a small scaling range still observable. The nu- 
merical investigations on the Cayley tree are described 
in detail in section 3. 

We measured directly the exponents r||, 7 and D. 
These exponents are presented in table 1. Since only 2 
exponents are independent in each case, measuring three 
different exponents provides a consistency check. The 
values displayed in table 1 do obey the scaling relations 
(y) and (|y), within numerical error. 

The most remarkable feature in table 1 is the mono- 
tonic decrease of Xc = -D — d as a function of d. In 6 
dimensions the APs become flat to within numerical ac- 
curacy, indicating that dc '^ &■ This is in accord with our 
previous estimate for the value of dc based on dynamic 
scaling. Since it has been also confirmed in |T^ that 
z = Zper we believe that our conjecture about the nature 
of the dynamical process of invasion is very plausible. We 
also note that in 4 dimensions the measured Xc becomes 
smaller than Xa, the annealed roughness. This implies 
(cf. eq. (S)) that the interface roughness x is dominated 
by annealed mechanism. In the next subsection we dis- 
cuss this point and analyze in detail its consequences. 

Another feature evident from table 1 is that in 6 dimen- 
sions the scaling exponents seem well converged toward 
their infinite dimensional values. 



TABLE I. Results of numerical simulation of bounded 
slope models. The finite d entries refer to simulations on hy- 
per-cubic lattices, and the d = 00 entry refers to simulations 
on a Cayley tree with coordination number 3. The exponents 
D, T were directly measured, as well as 7 in the finite di- 
mensional case. The value of r was calculated from D and ry 
using (W), in the finite dimensional simulations, and measured 
directly in the Cayley tree simulation. The value of x (the 
overall roughness) was not measured but is given for refer- 
ence. The numbers in parenthesis are errors in the last digit. 
We have checked that the numerical values of exponents pre- 
sented here satisfy all the theoretical scaling relations which 
connect them. 
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5.05(5) 


1.68 


0.2 


00 


.150(1) 
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6. Overall interface roughness 

We now turn our attention to the scaling of the growing 
interface as a whole. Above dc the APs are flat, and the 
invasion process becomes very similar to the annealed 
process: Each AP activates each member of a set of sites 
once, or at most a few times. In this regime the interface 
may be described by the annealed KPZ 



dth{x, t) = V2/i(x, t) + A(V/i(x, t))2 + 7/(x, t), (14) 

where h is the height and the random noise r/ describes 
the activation of sites by the flat APs. Since the acti- 
vation of sites by APs favors nearby sites, rj has spatial 
correlations derived from the distribution and structure 
of the flat APs. We expect that in large enough dimen- 
sions, the spatial correlations of 77 will become irrelevant, 
and that the critical exponents of the interface will iden- 
tify with those of KPZ with delta-correlated noise. We 
thus define the interface critical dimension di, such that 
Xa = Xq foi' d > di. When observing the invasion pro- 
cess of the interface as a whole, rather than the APs, 
di signifies the dimensionality where the effects of noise 
quenching become irrelevant. 

It is not difficult to estimate the statistics of the corre- 
lated noise 77, when d > dc, in terms of the AP properties. 
Since the APs are not correlated (for /o — fc, see 0), ?7 
is Gaussian, and its value at two points is correlated only 
if both sites were activated in the same AP. Therefore the 
correlation function may be estimated by the probability 
of such an event, 



(ry(x)?7(y)} - P(r|| > |y - x|)|y 



\D-d 



_ \D{2-T)-d 



(15) 



We now analyze the effects of those correlations, us- 
ing fc-space representation of the problem. We make the 
following definitions: Let 



h{'k,u)= f d'^xdth{x,t)e^''+"^*, 



(16) 



rjCkjUj) defined similarly. The response function is 

G{k, uj)6{k + k')5{uj + u') = {6h{k, uj)/6'n{k' ,uj')) (17) 
and we assume it has the scaling form 

G{k,cj)^k-^g{u;/k^), (18) 

where z is the dynamical exponent of annealed growth. 
The noise correlation in fc-space is 

(77(k, u;)Tj{k', J)) = fc"(5(k + k')5{Lo + Lo'), (19) 

where a = D{t — 2). 

The roughness of an interface described by eq. (Ilj) 
is generated by a combination of the noise correlations 



and the intrinsic roughness stemming from the nonlinear 
process. It is possible to give a naive calculation for the 
correlation induced roughness which is a lower bound to 
the actual roughness. To this end we define the "bare" 
height field as 

hoiKuj) = G{k,uj)T^{k,uj) (20) 

The bare correlation function is 

{ho{k,oj)ho{k' ,uj'))^G{k,Lo)G{k' ,Lo') {r^{k,uj)i^{k' ,Lo')) 
- k~^'-+°'g{Lo/k''f5{k + k')5{uj + uj'), 

(21) 
and the simultaneous bare correlation is given by 



_ 7,-2z + Q 



(/io(k,t)/io(k',i)) = fc 



so we get the bare roughness 



dujg(uj/k''Y ^ k 



-z+a 



2xo — max(z — a — d,0). 



(22) 



(23) 



Thus, the correlation induced roughness dominates the 
intrinsic roughness unless 



2x > z — a — d. 



(24) 



One can check that the condition for irrelevance (gj) is 
always satisfied for uncorrelated noise (a = 0), and that 
it also reproduces the known condition for irrelevance of 
the noise correlation in 1 dimension, a > —1/2 ||l[. 

The condition for irrelevance (gj) was derived for d > 
dc, and is satisfied in this range, see table 1 and ||l^ for 
values of x and z for high dimensional KPZ; we can thus 
conclude that di < 6. The measured value of Xc becomes 
equal to numerical measurements of Xa at 4 dimensions 
within numerical accuracy, so di > 4. Since Xc <?C 1 in 5 
dimensions, it is probably safe to use the criterion (p3) 
also for d ~ 5. The criterion is satisfied so that we get 
the final estimate 



4 < d, < 5. 



(25) 



Independent measurements of the overall interface rough- 
ness [^ indicate that this is indeed the case, see table I. 
We conclude that above 4 dimensions the roughness ex- 
ponent of the overall quenched growth becomes identical 
with the annealed growth exponent. 



III. CAYLEY TREE 

In many statistical physics problems, the infinite di- 
mensional limit is obtained by defining the model on the 



Cayley tree. In problems with an upper critical dimen- 
sion the Cayley tree solution gives the mean field expo- 
nents. 

The recursive structure of the Cayley tree simplifies the 
problems because of lack of self interaction, and allows 
in many cases to obtain a closed form solution. A well 
known example is the solution of percolation on the Cay- 
ley tree [M. This solution is instructive since it reveals 
some of the characteristics of the percolation clusters. 
We are interested in the analogous problem for interface 
growth. In this case, however, it will be apparent that 
Cayley tree realizations of this model might reveals prop- 
erties which are sub-dominant for any finite dimensional 
realization of the model. We shall show that strong fluc- 
tuations may potentially destroy the usual kind of phase 
transition and may lead to anomalous critical properties. 
The Cayley tree realization which we define here displays 
a simple critical behavior with a strong background of 
anomalous behavior, but in other examples the anoma- 
lous behavior dominates. 

We define the invasion process on the Cayley tree as 
follows: The interface is defined by assigning a height to 
each site of the Cayley tree. As in the finite dimensional 
case, there is a quenched random force / which depends 
on the site in question and the height of the interface. A 
local slope is defined between every pair of neighboring 
sites, and the bounded slope condition prevents the local 
slope from increasing beyond a critical value by advanc- 
ing the sites which violate this condition. As usual there 
is a critical value of force /c, above which there is as finite 
probability to create an infinite cluster. 

The avalanche propagates upwards and outwards from 
the original site. Since we expect critical avalanches to be 
flat, propagation of avalanches will be mainly in the out- 
ward direction. However, two additional processes com- 
plicate the dynamics. First, the bounded slope rule im- 
plies that repeated activation of a single site is followed 
by activation of all its near neighbors, and the number 
of these neighbors grows exponentially with the number 
of repeated activations of the original site. Second, since 
the activations are isotropic, sites may be activated, cre- 
ate sub-avalanches, and then be reactivated by a sub- 
avalanche created by of their children. Such backward 
activations creates self-interaction in the dynamics, and 
is one of the main reasons for the difficulty in analyzing 
the model in a quantitative way. This self-interaction 
is the reason that the model on a l-dimensional chain, 
being a special case of the Cayley tree, is not trivially 
solvable. 

The general origin of the anomalous critical behavior 
in the size distribution of the clusters on the Cayley was 
discussed briefly in the previous letter |lQ] . We will give 
an example below, where this mechanism can be studied 
analytically, but flrst we would like to present the quali- 
tative arguments, since they have an essential role in the 
behavior of the phase transition. 



Consider an interface that is initially flat, and a growth 
process initiated from a single site with a fixed value of 
/q. There is a finite probability that the same site will 
be activated once more, activating the q nearest neigh- 
bors obeying the bounded rule, h subsequent activations 
of the same site will result in activation of the q{q — 1)'* 
neighboring sites, again as a result of the bounded slope 
condition. We observe that the number of activations in 
this sub-process grows exponentially with h. However, 
since the probability for each of the activation is /o, the 
probability for such an event is fg which is exponentially 
small. Thus, even though the probability to create such 
a height difference decreases exponentially with h, the ef- 
fect might become important because of the exponential 
avalanches that follow from such events. This effect can 
dominate the probability distribution of s. 

The actual number of sites which are activated by the 
original site of height h is larger than the previous esti- 
mate, since other sites can be activated by the exponen- 
tial sub-avalanche. Another complication arises since the 
actual interface is not necessarily ffat. However, since it 
is statistically flat, a sufficiently long sequence of acti- 
vations should have an exponential effect. Therefore we 
estimate the mean number of activations s{h) resulting 
from the event that a single site is of height h, by 



Kh) - /i^ 



(26) 



with /i > q—1. On the other hand, the probability that a 
site in an /g avalanche reaches a height h should decrease 
like /'* with / larger than, but of the order of /q. We 
can combine these two estimates to find the probability 
of an avalanche of size s. 



P[s{h) <s<s{h + l)]^ f" 



it follows that 



P{s) 



„(log//logM)-l 



rifo) 



(27) 



(28) 



Thus on the Cayley tree, the tail of the distribution of 
the APs sizes is always a power law, with a non-universal 
/o-dependent exponent. 

Near the phase transition there are two possible sce- 
narios. Either the transition large clusters are created by 
the present mechanism, or they are the result of critical 
ffuctuations. In the first case there is no typical scale for 
the mass of the clusters, as follows from (P8|). However 
there is a typical length scale which remains finite in the 
transition. 

The second case may be realized if for sufficiently small 
l/o ^ fell the probability to create a very large critical 
cluster is larger than the probability displayed in (pq). 
Thus, there will be a crossover scale, beyond which nor- 
mal critical behavior dominates. In particular, there will 
be typical scales for the mass and (chemical) length of the 
clusters, which diverge at /c, and the clusters will have 



a finite fractal dimension. Since the mechanism which 
creates anomalous clusters is not universal, the question 
which critical behavior is realized depends on the specific 
realization. 

In the following subsections we are going to discuss, in 
addition to the model which we have already presented, 
two additional realizations of the bounded slope model 
on the Cayley tree. One realization was defined and an- 
alyzed in |pl|, and in it the interface grows into the tree, 
instead of in an additional direction. The other realiza- 
tion, to be defined in detail below, is a directed version 
of the model defined above, where backward activations 
are disabled. These two alternative realizations are in- 
tegrable, and display anomalous critical behavior. The 
original model, however, displays numerical evidence for 
crossover to standard critical behavior. 



1. Numerical analysis on the Cayley tree 

We performed Monte-Carlo simulations of the bounded 
slope model on the Cayley tree with coordination num- 
ber 3. We used the constant driving force algorithm ver- 
sion, i.e., for each site which was a candidate for activa- 
tion, a random number between and 1 was compared 
with a constant force /o, to decide whether the site should 
be activated. This process was repeated until either the 
cluster was blocked, or exceeded an upper bound in size. 

We collected the statistics of clusters, in terms the clus- 
ter size s, and in terms of maximal chemical distance from 
the source site i. In Fig. 1 we give examples of the size 
distribution for several values of /q. We observe three 
qualitative types of behavior, which can be explained in 
terms of the interplay between critical and anomalous 
fluctuations described above. In particular, the size dis- 
tribution very near to the critical value of forcing displays 
a crossover behavior, which we interpret as emergence of 
normal critical behavior. 
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FIG. 1. The distribution of cluster sizes n{s) multiplied 
by the critical distribution s'^" (arbitrary units), given as 
a function of s. We present distributions for forcing levels 
/o = 0.141, 0.145, 0.148, 0.150 ~pc and 0.155, displayed with 
higher values of / below lower ones. For low values of /, below 
Pc one can observe pure anomalous power law distributions, 
while as / increases a new feature becomes more and more 
apparent for large s, as one approaches Pc- The new power 
law seems to be independent of / near pc- Above Pc an expo- 
nential cutoff appears for large s. 

The distribution functions for subcritical driving force 
display a clear power law tail without a cutoff, but with 
an /o dependent exponent t(/o) (see eq. (p^)). The 
clusters in the tails of the subcritical distributions were 
created by the mechanism described in the last section. 
An unusual property which follows from the asymptotic 
form of the subcritical distributions, is the divergence of 
the average cluster size (s) when r < 2. This property 
has no significance beyond the fact that the size distri- 
butions have no cutoff in the subcritical range. 

The supercritical size distributions, on the other hand, 
display well defined cutoffs. The presence of these cutoffs 
stems from exactly the same reasons as the supercritical 
cutoff in usual critical phenomena; once a cluster reaches 
a certain size, it becomes very unlikely that it will not 
become an infinite cluster, and this is true regardless of 
the of mechanism which creates the power law tails in 
the subcritical range. 

Very close to pc (/o = 0.150) a third type of distri- 
bution is observed. There is sub-asymptotic power law 
tail with an exponent t(/o) characteristic of the subcrit- 
ical distributions, but at a certain value of s there is a 
crossover to a different power law distribution with an 
exponent Tc = 1.75 ± 0.05. We attribute the crossover 
to the emergence of normal critical behavior, which is 
masked by the anomalous fluctuations for other values 
of /q. The fact that very large clusters near /c are gen- 
erated mainly by critical fluctuations signifies that the 
transition itself is of normal type. 

Analysis of the length distributions of clusters p(^), dis- 
played in fig. 2, supports the picture obtained from the 
analysis of mass distribution. The subcritical distribu- 
tions decay faster than a power law, but for /o ^ /c the 
distribution decays asymptotically as a power law with 
a characteristic exponent ti — 3.8 ±0.1, consistent with 
normal critical behavior at fc- 



10" 



10' 



10 



10 



10° 



10" 



10' 



10 



10' 




1000 



FIG. 2. The distribution of clusters according to their 
chemical length i, multiplied by the critical distribution P'^ , 
for values of / corresponding to those displayed in fig. 1, 
except the smallest value. One can observe clearly that the 
distribution falls faster than a power for f < pc, and is power 
like at pc- 

We conclude that depinning transition of the bounded 
slope model on the Cayley tree is a critical phenomenon 
with a very strong background of anomalous fluctua- 
tions. Using the critical exponents Tc and ti one may 
also calculate the fractal dimension of critical the clus- 
ters £) = 2(t£ - 1)/(tc - 1) ~ 7. The large value of D 
shows that this phase transition is not in the percolation 
transition universality class. The critical exponents are 
also listed for reference in table 1. 



2. A solvable model 

It was remarked above that the main obstacle for ob- 
taining an exact solution for the invasion model on the 
Cayley tree is self-interaction through reactivation of a 
site by a "child" site. A natural way to overcome this 
difficulty is to define a model in which there is a preferred 
direction, 'back', in which the bounded slope rule is not 
enforced. The mechanism which creates anomalous large 
fluctuations is present in this simplifled directed model, 
and a more precise analysis of this feature is possible. 
However, in contrast with the "isotropic" Cayley tree 
model which was discussed in the previous subsection, we 
show that in the directed model anomalous fluctuations 
dominate up to the transition, and normal critical behav- 
ior is not observed. In addition to the directed model, 
we also discuss briefly the model presented in |1J], and 
reach similar conclusions. 

The directed model is therefore defined by precisely the 
same rules as the original model defined above, except 
that the slope may be arbitrarily large when the lower 
point is closer to the root. Thus defined, the process can 
be carried out also in the following manner: Pick a root 



site on the tree and let it be activated repeatedly un- 
til further activation is blocked by a pinning force larger 
than some fixed driving force pq. The probability that 
a height h will be reached by this process is Pq(1 — po)- 
Since the bounded slope rule works only forward, the 
interface height determined by this process may not be 
changed anymore, and is statistically independent of the 
interface height at any other site. Next, the same pro- 
cess is carried out for each of the child sites of the root, 
starting from height h — 1 ii h > 2; otherwise the in- 
vasion process stops. The process continues analogously 
for each of the children sites either until all the sites are 
blocked, or indefinitely, creating an infinite cluster (this 
can happen only when po > Pc)- 

The present model is equivalent to a branching pro- 
cess |15| ] , and we may use its recursive structure for anal- 
ysis, similar to the analysis of Cayley tree percolation. 
However, the model defines a universality class different 
from percolation, since it is a branching process with an 
infinite number of states: Each possible value of the inter- 
face height corresponds to a different branching process 
state. 

The mass distribution of clusters of the directed model 
may be calculated as follows. We denote the probability 
that the total mass of a cluster is s given that the inter- 
face height at the root site is h by p{s\h). This function 
also describes the probability that, given a child site with 
interface height /i, it will generate a sub-cluster of mass 
s. Since the total mass of the cluster is equal to the sum 
of the masses of all the sub-clusters plus the number of 
activations at the root, we have 

Pi^\h)= E E P{^i\hi)p{hi\h) 

si+...+Sg+h=s hi...hg (29) 

'Xp{s2\h2)p{h2\h) . ..p{Sq\hq)p[hq\h). 

The RHS of eq. ( p9| ) is a sum of products of expressions 
of the form 



f{s\h-l)^Y.P^s\h')p{h'\h) 

h' 

= Y,{l-p:,)p^^p{s\h' + h-\). 



(30) 



h'=0 



It is possible to invert (30) to express p(s\h) in terms of 
/, and rewrite eq. ( p9| ) 

(l-po)-^(-po/(s|/i+l) + /(s|/i)) 

E f{si\h-l)f{s2\h~l)...f{Sq\h~l). 



si+...+s„+h=s 



(31) 



The Laplace transform of the last equation is 
(1 - pQ)-\-pog{x\h + 1) + g{x\h)) =. x''g{x\h - 1)«, 

(32) 



where g{x\h) = J^s ^'^ fiM^) i^ the Laplace transform of 
f{s\h). The solution g of this non-linear difference equa- 
tion contains all the information about the mass distri- 
bution. In particular many properties of the distribution 
p{s\h) follow from the behavior oi g{x\h) near x = 1. For 
example, the probability to create an infinite cluster is 



Pooih) = 1 



g{l\h)-pog{l\h+l) 



(33) 



1 -Po 

Similarly, the expected cluster size starting at height h is 
d^g{l\h)-pod,g{l\h + l) 



{s\h) 



1 -Po 



(34) 



{s\h) obeys a linear difference equation obtained by lin- 
earizing (B3) near g{l\h). For small enough po, g(l|^) 
must be 1, and we can linearize eq. ( p2| ) near a: = 1 to 
get. 



-pod,gil\h+l) + d.,gil\h) 
(1-Po) 



qd^gil\h - 1) + h, (35) 



which can be solved in closed form, yielding expressions 
for {s\h). The solution ceases to be realizable when po = 
Pc, where Pc is defined by 



2qpc{l -pc) = 1- 



(36) 



This implies that for po > Pc 5(1 1 ^) ¥" Ij so that pc 
indeed is the critical value of po . 

Further properties may obtained by solving eq. (p3) 
numerically. One finds that for pq near pc and 1 — a; 
small 

g{x\h) - 1 + (1 - x)gi{h) + g,ih)il - x)^ , (37) 

where 2 < t < 3. This type of asymptotic behavior im- 
plies, using the Tauberian theorems ||ll|] |16|, that f{s\h) 
and therefore also p{s\h) behave asymptotically as s^'^ 
for large s, and in particular that {s\h) is finite but \s'^\hj 
diverges for pq < pc (near pc). For example, for g = 2, 
Po = Pc we find that r = 2.2 ± 0.05 

We observe that the subcritical mass distribution of 
the directed model is quite similar to that of the original 
model, with a po-dependent power law tail, but in con- 
trast with the original model, this behavior persists for 
values of Po up to and including pc- This indicates that 
the directed model belongs to the second case outlined in 
the beginning of this section, where the phase transition 
is dominated by the anomalous fluctuations caused by 
exponentially large compact avalanches. The same con- 
clusion also arises from the analysis of length distribution 
of the clusters, presented next. 

The analysis of cluster length is more involved than 
that of cluster mass. Using the analogy with an infinite 
state branching process, we interpret the sites at a given 
chemical distance i from the root, as the £th generation 



of the process. The population of the process is defined 
by the values of the interface height, where each height 
corresponds to a different species of the branching pro- 
cess. Thus, we may define the ^th generation population 
vector Zi — {Z\ , Z] ,....), where Z\ ' gives the num- 
ber of sites at chemical distance I from the root with 
interface height equal to 1, etc. The branching process is 
defined by giving the probability p„(Z) that an individ- 
ual of type n will generate a set Z = {Z^^\ Z'^'^\ ....) of 
individuals in the next generation. The values of these 
probabilities follow from the definition of the directed 
model. Since we were mainly interested in the asymp- 
totics of the length distribution of clusters, we analyzed 

the dependence of the moments Izf^j and iZ^Zg^ 
on t. Such moments average both over clusters shorter 
than (., contributing 0, and over clusters longer or equal 
to i. For critical clusters we expect that the distribu- 
tion of Z^ does not depend on the cluster size, given that 
cluster is longer or equal to (.. Thus, we may write 



Z, 



(") 



,(n 



= p(£)(Z)"^|Z, ^0)^c„p(£)s(^) 



(38) 



where p{i) gives the probability that a cluster is longer 
than £, and s(^) is the average number of sites at a dis- 
tance i from the root. s(£) would be of the form lf^~'^ 
if the critical clusters have a finite fractal dimension Dg. 
Similarly 



ry{n) y{m) 



n.p{i)s{lf 



(39) 



so that we may deduce the behaviors of p{€) and s{t) 
from those of the first and second moments of Z. 

Explicit calculations, not reproduced here, show that 
first moments behave as a power law for large ^, while 
second moments grow exponentially. This implies, using 
(p8| ) and (p9|), that s(t) grows exponentially and p{i) de- 
cays exponentially with L In other words, the clusters 
have a characteristic length, and have an infinite fractal 
dimension. These properties are expected for a dcpinning 
transition dominated by anomalous fluctuations, as was 
already indicated by the analysis of cluster mass distri- 
bution. We conclude that in the directed model critical 
fluctuation are too weak to be observed even near pc- 

An alternative version for defining the bounded slope 
model on the Cayley tree was given by pT| . In this ver- 
sion, the Cayley tree is divided into sets of sites with 
a fixed height, and the invasion process occurs directly 
from one site to its neighbor, rather than by raising the 
interface in a single site as in the version presented in 
the present paper. The paper iTlt] presents an analysis 
of the cluster mass distribution. As in the models pre- 
sented here, the mass distribution displays non-universal 
power law asymptotics, stemming from anomalous fluc- 
tuations. Since there is no indication of a crossover to a 
different type of asymptotics near pc^ we conclude that 
in this model the transition is anomalous. 
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3. Consequences for large but finite d 

In many cases one uses information from Cayley tree 
realizations of statistical models to deduce properties of 
high-dimensional realizations. In particular, one may 
hope to find "mean field" behavior which should be valid 
above a finite critical dimension dc, if such exists. 

In the case of the bounded slope models, however, we 
have a dilemma: The behavior on the Cayley tree is com- 
plicated by the creation of exponentially large avalanches 
so we observe two types of critical behavior. In some 
cases we had a fiuctuation dominated phase transition 
and in others the phase transition is of standard critical 
type. 

We believe that the properties observed on the sym- 
metric Cayley tree should serve as a better guide to the 
finite dimensional case, for two reasons. First, the sym- 
metric model we used is closer in its definition to high di- 
mensional bounded slope models than both the directed 
model we solved and the tree model suggested in [ pj , 
since it has a definite height and symmetric activations. 

The second reason follows by examining the process 
which creates the anomalous avalanches in large but fi- 
nite d. The analogous arguments which led to the con- 
clusion that the cluster mass distribution has power law 
asymptotics (cf. eq. p8[ ) repeated in the case of finite 
d yield a distribution which behaves asymptotically as a 
stretched exponential exp(— s^/'*). Such behavior will al- 
ways become sub-dominant for pq very close to Pc- Thus 
in a finite dimension we always expect a crossover from 
an anomalous fiuctuation regime, to a critical regime, 
such as displayed by the symmetric Cayley tree model. 
We thus conclude that in large dimensions the associ- 
ated processes become objects with a fractal dimension 
approaching Dqo = 7 ± 1. Since in six dimensions D, as 
well as other critical exponents are close to their infinite 
dimensional values, we are led to conjecture that there 
exists a finite critical dimension beyond which infinite- 
dimensional critical behavior is achieved. 



the exponent Xg for d > A using this understanding. A 
basic conclusion of this is that different driving methods 
of this model will display different statistical properties. 

Another basic idea is that the dynamic process of 
avalanche formation below 6 dimensions is a percolative 
process. There are two main conclusions from this claim. 
Firstly, the roughness exponent Xc goes down with in- 
creasing dimension, and the dimension where Xc = is 
dc — 6. Secondly, the dynamic exponents of the process 
should be the same as in percolation [Q. 

From our analysis on the Cayley tree we find that in 
high dimensions strong fluctuations appear as a result of 
exponentially rare events, as well as a result of critical 
behavior. These fluctuations are observed by the pres- 
ence of wide tails in the AP size distribution, and by 
power law tails in the cluster distribution on the Cayley 
tree. As a result, one can observe two types of scenarios 
for the depinning transition on the Cayley tree. In ad- 
dition to a critical transition, it is also possible to have 
a fiuctuation dominated transition, characterized by a 
correlation length which does not diverge. Indeed, this 
behavior is observed both in a directed model we solve 
and in a model analyzed in |ll| . On the hand, on a 
simple symmetric bounded slope model the usual critical 
transition is realized. We conclude that the latter model 
is the appropriate high dimensional limit. 

Examination of the values of scaling exponents ob- 
tained numerically (see table 1) indicates that there is 
a finite absolute critical dimension, in which all the ex- 
ponents reach their infinite dimensional limit. 

In summary, we find that there are at least three qual- 
itative transitions in the basic structure of an invading 
interface when dimensionality is increased. This refiects 
the wealth of processes occuring during such an invasion 
process. 



IV. CONCLUSIONS 

The analysis presented in this paper allows us to draw 
conclusions on the processes belonging to the quenched 
KPZ universality class and processes of percolation of di- 
rected surfaces. We believe that studying the behavior in 
different dimensionalities has provided a lot of interest- 
ing insights. The main conclusions and ideas presented 
in this paper in our view are as follows: 

There exist two different roughness exponents charac- 
terizing the invading front: The exponent Xc of the AP 
roughness, and the overall interface roughness Xq- The 
overall roughness is created by random deposition of ob- 
jects with a roughness Xc- We are in fact able to predict 
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